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Abstract 

The two-orbital superconducting state is modeled by on-site intra-orbital negative- U Hubbard correlations together 
with inter-orbital pair-transfer interactions. The critical temperature is mainly governed by intra-orbital attractive 
interactions and it can pass through an additional maximum as a function of band filling. For the certain number of 
electrons the clear interband proximity effect is observable in the superconducting state of the band with a smaller 
gap. The influence of band fillings and orbital site energies on the temperature dependencies of two-component 
superconductivity coherence lengths and magnetic field penetration depth is analyzed. The presence of proximity 
effect is probably reflected in the relative temperature behaviour of characteristic lengths. 

Keywords: two-orbital superconductor, coherence lengths, magnetic field penetration depth, Hubbard model, 
interband proximity effect 



1. Introduction 

The multi-component microscopic nature of super- 
conductivity has been established in a number of ma- 
terials (MgB2 HI, cuprates Jptl, strontium rutenates 
(3], iron-arsenic compounds (UNb-doped SrTi03 @1> 
NbSe 2 El, V 3 Si J3l etc.). The first relevant theoretical 
scenarios have been discussed in the literature long time 
ago JESEl and later on they have been developed for 



Bill 111 111 iMlllJ and 



various systems (see e.g. Refs. 
references therein). 

The peculiar properties of multi-component super- 
conductors appear due to the involvement of coupled 
electron subsystems (orbitals, different Fermi surface 
sheets or other elements of electron structure) into the 
formation of superconducting ordering varying substan- 
tially the physical nature of the phenomenon and open- 
ing the possibilities for novel effects. As a result of the 
presence of inter-component pairing the certain quanti- 
ties (e.g. coherence lengths and the relaxation times of 
superconductivity fluctuations) cannot be attributed to 
the initially independent subsystems and they describe 
the joint features of the whole two-component conden- 
sate CIE H Cl- 
in a two-orbital superconductor with inter-orbital 
pairing, similarly to the results of Refs. |17| 20l 21 \, 
one of two solutions for coherence lengths diverges 



near the phase transition point as a function of tem- 
perature, while the other one is non-critical 111 811 . We 
note also that the existence of critical and non-critical 
length scales in the coherency associates with the con- 
ception of critical and non-critical fluctuations |[l9l I22I1 
in a two-component superconductor damping in differ- 
ent time-scales. 

In the present contribution we continue to examine 
the temperature behaviour of characteristic lengths of 
a two-orbital superconductor described by negative- U 
Hubbard model 02311 with inter-orbital pair-transfer in- 
teraction in dependence on the number of electrons. Our 
preliminary study has been done focusing on band fill- 
ing effects in the spatial coherency including Van Hove 
singularity 111 811 . However, in that work 111811 the or- 
bital energy differences and charge fluctuations were 
neglected. In the present paper we clarify their role in 
the scaling of critical and non-critical coherence lengths 
and magnetic field penetration depth by band fillings. 
The numerical calculations have been carried out for a 
two-dimensional square lattice. 



2. Superconductivity in a two-orbital system 



Two-orbital superconductivity is modeled by the 



Hamiltonian of the following form (see e.g. II24I0 : 
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where a + (a,-^) is the electron creation (destruction) 
operator in the orbital a — 1,2 localized at the site i; 
cr is the spin index; ff? is the hopping integral; e° is 
the orbital energy; // is the chemical potential; U aa < 
is the intra-orbital attraction energy; n m(T = at a ia(r 
is the particle number operator; U aa ' with a + a' is 
the inter-orbital interaction energy. Note that the both 
intra- and interorbital interaction channels involved lead 
to superconducting ordering. 

In the reciprocal space the Hamiltonian (03 can be 
represented as 

H = [ £ff(k) " ^ a «kcr<W 

a k tr 

+ 2JV Z Z Z Z ^ a A<T a t(-k+q)-<T 
a,a' k,k' q cr 
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where ^(k) is the electron band energy associated with 
the orbital a and N is the number of lattice sites (num- 
ber of atoms). The bulk superconducting state is de- 
scribed on the basis of the Hamiltonian © by means of 
the Hartree-Fock-Gorkov self-consistent equations 

A„ = N- l J]u aa 'J](a ^ kl a a ^) 
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Here E a (k) = ^(k) + |A„| 2 , s a (k) = e„(k) + 
5 U aa n a — fi, n a are the occupation numbers of orbitals 



(bands) per lattice site, and « is the total number of elec- 
trons per site. The Eqs. Q determine the superconduct- 
ing gaps of orbitals (bands), the Eq. © gives us the 
orbital occupation numbers and the Eq. (0, the chem- 
ical potential fi in dependence on n and other relevant 
quantities. The temperature of superconducting phase 
transition T c satisfies the equation 



l + U u gl (T c ) U l2 g 2 (T c ) 
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In the general case of spatial inhomogeneity and non- 
zero magnetic field we have the following Ginzburg- 
Landau equations 



A„(r) = -YjU 
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q=0 



A is the vector potential, j s is the density of supercur- 
rent, Oq is the magnetic flux quantum and Vq is the vol- 
ume of unit cell. According to the Eq. (fTTT i we restrict 
ourselves to the isotropic electron energy spectrum. 
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Figure 1 : The critical temperature T c versus band filling n = n \ + n% 
for various e?/f = 0, 1, 2, 3, curves (a)-(d), respectively. 



3. Characteristic lengths of a two-orbital supercon- 
ductor 

By considering the solutions of Eqs. © which cor- 
respond to the small deviations of gaps from their bulk 
values, r] a (r) = A„(r) - A a , one can find two length 
scales characterizing the spatial behaviour of supercon- 
ductivity in a two-orbital system: 

G(T)± VG 2 (T) -4K(T)y 
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The soft or critical coherence length % S (T) diverges at 
the phase transition point T = T c , while the rigid or 
non-critical coherence length g r (T) remains finite. 

On the basis of Eq. (|9]i one obtains in a standard way 
the magnetic field penetration depth for a two-orbital 
superconductor, 
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where jjQ is the magnetic permeability of free space. 



4. Numerical results 

Below we present the results of numerical calcula- 
tions for two-dimensional square lattice with hopping 
integrals between nearest neighbours ff" - t and elec- 
tron band energies associated with s-orbitals s a (k) = 
e° - 2t \cos(ak x ) + cos(afc v )J, a is the lattice constant, 
and — - < L, < -. The intra- and interorbital in- 
teractions are chosen as U u = -1.5/, U — -2.5f, 
\U l2 \ = \U 2l \ = 0.04?, and the site energy e 2 is fixed 
at e° = 0. In all calculations we have chosen ks — 1 . 

Using Eqs. © and (0 we calculated the supercon- 
ducting critical temperature T c . Fig. 1 shows the de- 
pendence of T c on the band filling n = n\ + n%. One 
can see that the curve T c vs n is most symmetric as 
£j = e° ( see Figl> curve (a)). In this case, the small 
shift of the maximum of T c towards the smaller values 
of occupation number n is caused by effective Hartree 
corrections U aa n a /2. First of all, these corrections are 
different for various orbitals as the intra-orbital interac- 
tions U aa are different. Consequently, the distances be- 
tween the chemical potential fi position and the lower 
edges of orbital bands are slightly different. There- 
fore the position of the maximum of T c is moved re- 
flecting the above mentioned difference which appears 
here due to the choice of intra-orbital pairing potentials 
\U 22 \ > \U n \. 

By increasing the difference between the bare site en- 
ergies |gj - e°l we observe that the asymmetry of the 
function T c (ri) increases. In the cases (a) and (b) in Fig. 
1, the site energies of orbitals are closer and the influ- 
ence of the Hartree terms does not lead to dramatically 
different distances between the chemical potential and 
lower band edges. As a result the single peak in T c vs n 
is preserved. However, with the further increase of the 
difference between of bare site energies 9 the effect of 
(15) the Hartree renormalization \e®+U 11 n\ l2—e^—U 22 ri2l2\ 
becomes stronger and we observe the splitting of the 
single peak of T c into two maxima (see Figl, curves (c) 
and (d)) caused by the circumstance that the chemical 
potential passes in the region < n < 4 two Van Hove 
singularities related to the bands a — 1,2. It should 
be also noted that the larger values of - e^l mean 
the larger difference between the band occupation num- 
bers «i and «2 which corresponds to the stronger Hartree 
renormalization. 

Next we solve the gap equations (Eq. Q) below T c . 
In Fig. 2 one can see the plotted superconducting gaps 
A a vs temperature for different band fillings. Note that 
the increase of n changes the order of Ai and A2 in the 
energy scale for fixed temperature, c.f. Figs. 2a, 2b 
and 2c, 2d. With that for the intermediate values on n 
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Figure 2: 

fixed £j/f = 2 and n 
tively. 



The temperature dependence of superconducting gaps for 
1.10,2.57,2.64,2.90, panels (a)-(d), respec- 



the gaps Ai and A2 are relatively close (Figs. 2b and 
2c) while for n = 1.1 and n = 2.9 one of the gaps 
strongly dominates (Figs. 2a and 2d). The latter pecu- 
liarity is caused by the chemical potential vicinity to the 
Van Hove singularity which supports superconductivity 
in the corresponding band. The weaker superconduc- 
tivity in another band is related, at least partially, to the 
interband proximity effect. At the same time the band 
fillings n = 1.1 and n = 2.9 correspond approximately 
to the maxima of T c (n) in Fig. 1. Thus, the formation 
of two peaks of T c as a function of n in Fig. 1 reflects 
the redistribution of the driving roles of bands in the ap- 
pearance of superconductivity. 

In a single-orbital (single-band) system the Ginzburg- 
Landau parameter k = A/g defines the type-I (k < 
1 / y2) and type-II (k > 1 / V2) superconductivity. The 
coherence lengths (see Eq. ([T3T l) and magnetic field pen- 
etration depth (see Eq. ( TTTl i) for two-orbital model are 
plotted in Fig. 3 for the band fillings discussed above. 
Apart from the divergence of the soft length scale tj s 
at T c , we observe also a maximum for the rigid length 
scale £ r slightly below T c in Figs. 3b and 3c. This is due 
to the closeness of the superconductivity gaps for these 
band fillings, see Figs. 2b and 2c. 

Another interesting observation is related to the 
crossing point of and A. In Figs. 3b-d this cross- 
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Figure 3: The temperature dependence of coherence lengths and 
magnetic field penetration depth in units of the lattice space a for fixed 
e\lt = 2 and n - 1. 10, 2.57, 2.64, 2.90, panels (a)-(d), respectively. 



ing appears at the temperatures very close to T c while 
in Fig. 3a the crossing point is located remarkably 
lower in the temperature scale. Thus, there exists al- 
ways a small temperature region in Figs 3a-d where 
£ r < A < f s Q. One can treat the fulfilment of the in- 
equality £ r < y[2A < as necessary but not sufficient 



We have used the relatively small values of attractive interaction 
\U aa \ to make the Hartree-Fock-Gorkov approximation plausible. In 
this situation we found that the penetration depth in our calculations 
is always smaller than the critical coherence length. The negative-U 
Hubbard model for larger \U aa \ has to be treated in a different way 
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condition for type- 1.5 superconductivity 026II . see Ref. 
lfl6Tl . To shed more light on the problem one has to study 
the character of the interaction between vortices [16]. 

The larger extent of the temperature region in Fig. 
3a where % r < A < f s may be related to the interband 
proximity effect which is more pronounced in this case 
(see Fig. 2a). It was shown in 12711 that type- 1.5 super- 
conductivity can arise due to the proximity effect in a 
two-band system. 

5. Conclusions 

We presented the results of band filling effect on two- 
orbital superconductor. The dependencies of critical 
and non-critical coherence lengths on temperature were 
compared to the penetration depth for different electron 
numbers. It was shown that the band structure including 
site energy and position of Van Hove singularity is im- 
portant for superconducting gaps temperature behaviour 
as well as in the relationship between the magnitudes 
of critical and non-critical coherence lengths and pene- 
tration depth. We have identified a temperature region 
where . < A < f s . Its larger extent could be associ- 
ated with the presence of interband proximity effect in a 
two-orbital superconducting system. 
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